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Abstract. We performed a spectral analysis of the EUVE 
spectrum and ROSAT data of the unique object H1504+65, 
which is an extremely hot post-AGB star entering the white 
dwarf cooling sequence. It is the only pre-white dwarf known, 
whose surface is free of hydrogen and helium, hence, it repre- 
sents the bare core of a former AGB star. The EUV spectrum 
(75-150A) is dominated by strong O vi lines and we can iden- 
tify a number of Ne vn lines. 

EUVE and ROSAT data can be fitted with models between 
T Gff =170000K and 200 000 K. We derive an extraordinarily 
high neon abundance (2%— 5% by mass, i.e. 20-50 times so- 
lar) which we confirm by an optical Ne vn line detected in 
a Keck HIRES echelle spectrum. This abundance is expected 
for 3a processed matter and corroborates our understanding of 
HI 504+65 as a C-O stellar core which has lost its entire H- and 
He-rich envelopes. 

Key words: Stars: abundances - Stars: atmospheres - Stars: 
evolution - Stars: AGB and post-AGB - white dwarfs - Stars: 
individual: H1504+65 



1. Introduction 

The optical counterpart of the soft X-ray source H1504+65 was 
detected by Nousek etal. (1986) and spectroscopy revealed 
the unique nature of this object. From the spectroscopic sig- 
natures it belongs to the PG1159 stars, which are very hot 
(T off =75 000-1 80 000 K) hydrogen-deficient pre-white dwarfs 
(log 5=5.5-8, cgs-units). But in contrast to all other PG 1159 
stars it appeared to be also helium-deficient, which was later 
confirmed by model atmosphere analyses of optical spectra 
(Werner 1991) and FUV spectroscopy with the Hopkins Ul- 
traviolet Telescope (Kruk & Werner 1998). Detailed NLTE 
line profile fitting revealed that HI 504+65 is the most mas- 
sive PG 1159 star, having the highest surface gravity, and that 
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it is among the hottest of this group (M/Mq = 0.86 ± 0.15, 
log 5 =8.0± 0.5, T cff =170000K± 20 000 K; Werner 1991). 

HI 504+65 is one of the brightest sources in the EUV and an 
early attempt to understand its spectrum recorded with the Ex- 
treme Ultraviolet Explorer (EUVE) failed in several respects, 
particularly the model flux was overestimated by an order of 
magnitude, however, it became obvious that the spectral ap- 
pearance is dominated by strong lines of O vi (Barstow etal. 
1995). Since that time considerable progress was made in 
NLTE modeling of stellar atmospheres so that a new attempt 
to analyze these data together with soft X-ray ROSAT observa- 
tions and new optical Keck spectra seemed promising. 

2. Observations and data reduction 

The EUVE spectrum of HI 504+65 was retrieved from the 
public archive. The exposure lasted 38 000 seconds, starting 
on Dec. 5, 1993. It was reduced (including flux calibration 
and subtraction of higher order contributions) with the stan- 
dard procedures of the IRAF/EUV software package (Version 
1.6.2). H1504+65 is only detected in the EUVE short wave- 
length spectrometer because of its relatively high interstellar 
column density. The observed spectrum covers the range 75- 
150A at a resolution of 0.5 A. 

A useful complement are X-ray ROSAT PSPC data, which 
were retrieved from the public archive. The observation lasted 
4882 seconds, beginning on July 18, 1990. A pulse height dis- 
tribution was extracted by using the MIDAS-EXSAS software. 
Source flux is detected in the range 30-120A (100^120eV). 

We have performed optical spectroscopy of HI 504+65 at 
the 10m Keck I telescope on July 20, 1998, under excellent 
weather conditions. We have taken two spectra, which were 
later co-added, with a total of 6000 seconds exposure time, us- 
ing the high resolution echelle spectrograph (HIRES) and the 
blue cross disperser. For details on data reduction see Zucker- 
man & Reid (1998). The spectrum covers the full wavelength 
region between 3600A and 5100A at a resolution of 35 000. 

3. Model atmospheres and spectral fitting 

We have computed a grid of non-LTE model atmospheres with 
different input parameters (T c g, \ogg, element abundances). 
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Fig. 1. Top panel: Emergent flux of a model with T o ff=175 000K, log g =8. Shifted upward is the same model spectrum but de- 
graded with a 0.5 A FWHM Gaussian. Overplotted (dashed line) is the continuum flux computed without occupation probabilities 
in order to display the position of the absorption edges. They are smeared out by atomic level dissolution and broad absorption 
lines. The spectrum is dominated by strong and broad O vi lines and less prominent O v and Ne vn lines. Bottom panel: Com- 
parison of the observed EUVE spectrum with this model (attenuated by an interstellar column density N(H i) = 5.1 • 10 19 cm -2 
with He n/H i = 0.068 and He i/H i = 0.052) 



The models include the most abundant elements, C and O, as 
well as Ne and He self-consistently and they are plane-parallel 
and in radiative and hydrostatic equilibrium. The computer 
code is based on the Accelerated Lambda Iteration method and 
is described in detail by Werner & Dreizler (1999). 



Let us present the atomic input data at some length, be- 
cause this is the first detailed analysis of an EUV spectrum of a 
PG 1 159 star with data of this good quality. We note in passing 
that similar models were applied only in the case of one other 
PG 1 159 star (PG 1520+525; Werner et al. 1996a) from which 
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a useful spectrum could be taken by EUVE. However, the S/N 
was much worse in that case. 



3.1. Potential spectral lines and absorption edges in the 
EUVE spectrum ofH1504+65 

No features of He n and C iv are located in the 75-150A spec- 
tral range, however, C iv absorption edges at longer wave- 
lengths are important because they contribute strongly to the 
background opacity in the EUVE region. Although C v is the 
dominant ionization stage of carbon we do not expect to see 
any of its spectral lines because temperatures are too low to 
populate excited C v levels. On the other hand, we know from 
previous studies that the C v ground state edge (3 1 .6A) almost 
completely blocks the flux (Werner etal. 1996b) in a PG 1159 
stellar atmosphere. This has no observable influence for the 
EUVE range, but is important in order to interpret correctly 
the ROSAT PSPC data. 

The strongest opacity source is oxygen. O vi is most impor- 
tant: In the spectral range in question we find two absorption 
line series arising from the ground state and the first excited 
state. The respective absorption edges are located at 89. 8A and 
98. 3A (see Fig.|J top panel), the latter one causes the flux to 
drop at A < 100A. Absorption edges of O v are weaker (be- 
cause O vi is more populated), but they are more numerous so 
that their combined opacity is important for the overall flux dis- 
tribution (Fig. [I]). These edges arise from the six lowest O v lev- 
els, where we have accounted for the possibility that from any 
2s2p configuration either the 2p or 2s electron can be ionized, 
leaving behind an O vi ion in the ground state or first excited 
state, respectively; therefore two edges appear at different lo- 
cations in the spectrum. Spectral lines of O v are less important 
because of the weak population of this ion and because the line 
profiles are intrinsically less broad than those of the O vi lines 
(quadratic vs. linear Stark effect, see below). 

Many lines from Ne vn are located in the 80-140A range, 
as well as numerous absorption edges from Ne iv to Ne vn. The 
edges turned out to be very weak as compared to oxygen, which 
is of course a consequence of the lower Ne abundance. 

Primary source for the used level energies is Bashkin & 
Stoner (1975). Carbon and oxygen oscillator strengths for spec- 
tral lines as well as bound-free cross-sections for photon and 
electron collisional ionization were largely provided by K. But- 
ler (priv. comm.). Oscillator strengths for Ne vn lines were re- 
trieved from the Opacity Project data base (Seaton et al. 1994), 
and bound-free cross-sections were computed hydrogen-like 
(the Ne model atom is essentially identical to that of Werner & 
Rauch 1994). Alternatively we have also used Opacity Project 
photon ionization cross-sections for all species, but the conse- 
quences are not relevant for our present study. 

Finally we note that helium, for which an upper limit of 1 % 
by mass was derived, does not at all affect the EUV spectrum. 



3.2. Line broadening and pressure ionization 

Line broadening is a problem which can be tackled only ap- 
proximately at best, particularly in the most important case, the 
O vi lines. O vi is a hydrogenic ion. Its energy levels with equal 
principal quantum number are closely spaced but not degener- 
ate. As a consequence, line broadening ranges between the lin- 
ear and quadratic Stark regimes. We faced the same problem for 
optical lines in previous analyses and we refer to Werner et al. 
(1991) for details on the approximative approach which we ap- 
ply here, too. For O v and Ne vn lines we assume quadratic 
Stark broadening, which is reasonable except for the highest 
members of the O v line series in our synthetic spectra, because 
linear Stark effects could become important here, too. 

Another problem is posed by the numerical treatment of 
the bound-free photon cross-sections. It is clear from the out- 
set that absorption edges in a stellar spectrum are not really 
sharp. Pressure effects lower the atomic ionization potential 
and tend to smear them out (see Fig.pl top panel). We account 
for this effect by an occupation probability formalism (Hum- 
mer & Mihalas 1988) which was generalized to NLTE con- 
ditions by Hubeny etal. (1994). We use essentially their nu- 
merical treatment, but two details demand special attention in 
our case. First, the perturber particles, which impose an electric 
microfield at the location of the radiating atom, are not protons 
(as it is the case in atmospheres with solar composition), but 
highly charged particles (mostly C v and O vn). Consequently 
the critical field strength f3 c as given in Hubeny etal. (1994; 
their eq. A. 2) needs to be scaled by C 1 / 3 , where C is the ratio 
of the total ion density to the electron density, n\ on /n c . And 
second, the plasma correlation parameter a, which is defined 
as the ratio of the mean distance of ions to the Debye length, 
needs to be scaled, so that we have instead of their eq. A.4: 
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where N\ and Z x are the number density and charge of the i- 
th ion, respectively, and n c and T are the electron density and 
temperature. This procedure is formulated strictly for hydro- 
genic ions only and, hence, applicable for O vi, but we also 
treat the O v edges in this way, because we lack any better ap- 
proach. But since the O v edges are weak, this uncertainty is 
not important for our analysis. 

To summarize, the most severe uncertainty in the synthetic 
spectrum calculations is the lack of a satisfactory treatment of 
O vi line broadening, because strong O vi line merging has a 
significant effect on the overall shape of the EUV spectrum. 
Theoretical line broadening data exist for only few of the O vi 
lines relevant here, but these do not account for linear Stark 
effects (Dimitrijevic & Sahal-Brechot 1992). 

3.3. Iron group elements 

EUVE spectra of hot hydrogen-rich white dwarfs are com- 
pletely dominated by iron and nickel lines and bound-free con- 
tinua (e.g. Wolff etal. 1998), hence, it was felt mandatory to 
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check the relevance of these opacities in the case of HI 504+65. 
Model calculations show that effects of the iron group elements 
(with a solar abundance fraction) are not detectable at the res- 
olution and S/N level of the present EUVE spectrum. Contin- 
uous opacities of the Fe group are negligibly small compared 
to the dominant species. On the other hand, our model predicts 
that a large number of lines should be detectable in future ob- 
servations with the Chandra X-ray observatory, which will be 
able to provide spectra with higher resolution and better S/N. 

3.4. Treatment of interstellar absorption in synthetic spectra 

EUV and soft X-ray spectra are very sensitive to absorption 
from interstellar hydrogen and helium. For the EUVE spec- 
trum this attenuation was calculated according to the model of 
Rumph et al. (1994). The interstellar column density of neutral 
hydrogen (iV(H i)) enters as a free parameter. It was chosen 
so that the observed flux could be reproduced at A ^ 130 A. 
The column densities of He i and He n were fixed relative to 
hydrogen using the mean values of He i/H i = 0.068 and 
He n/H i = 0.052 from Wolff etal. (1999). 

For the ROSAT PSPC pulse height distribution we have 
also used fixed relative abundances for helium (see Jordan 
etal. (1994) for details of the analysis). Since different spec- 
tral regions are used to determine the interstellar absorption the 
N(H i) values from the ROSAT and EUVE analyses differ. 

4. Results and discussion 

In the top panel of Fig. [y we display a model spectrum for 
HI 504+65. The numerous O and Ne lines are identified and 
close inspection shows that the O vi lines are strongest. Line 
blending with O v and Ne vn occurs at many locations and the 
degraded spectrum (shifted upward in the Figure) which simu- 
lates the EUVE resolution of 0.5A demonstrates that disentan- 
gling of single lines becomes difficult. However, it is possible 
to identify a few isolated lines or line cores of O vi and Ne vn 
which are indicated by arrows in the lower panel of Fig. [I]. The 
roll-over of the flux towards shorter wavelengths is caused by 
the strong absorption edge of the first excited level of O vi at 
98A, which can be seen by the dashed line in the top panel 
of Fig. [j] which represents the continuum flux of the model. 
As already mentioned, the sharp absorption edge is strongly 
smoothed out by pressure effects as well as the converging ab- 
sorption line series. For the same reason all other absorption 
edges which are indicated in Fig. fi] cannot be recognized in the 
final model spectrum. 

4.1. Effective temperature 

The best fitting model was found by taking models with dif- 
ferent T c ff and keeping fixed log g =8 and the abundance ratio 
C/0=1. The latter two parameters, when varied within the lim- 
its of the optical line analysis, affect the model spectrum much 
less than T e g, The model flux is normalized to the observed 
visual magnitude (V=16.24, Nousek etal. 1986) and the inter- 
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Fig. 2. EUVE spectrum as compared to two models with dif- 
ferent T cff : 175 000 K (dashed line) and 200 000 K (full line). 
The hot model fits the flux at A < 100A, but overestimates 
the maximum flux by a factor of two. Interstellar columns are 
7V(Hi) = 5.1 • 10 19 cm- 2 and 6.0 
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stellar column density, which is responsible for the roll-off at 
the long wavelength end of the spectrum, is treated as a free pa- 
rameter. The neon abundance can be estimated from the Ne vn 
line near 106A which is observed to have a deeper core than 
the blending O vi lines near 107 A. 

The best fit is displayed in the bottom panel of Fig. |[ with 
model parameters as indicated. The overall fit is satisfactory. 
The strongest O vi and Ne vn line features can be identified and 
the model can fit them. The absolute flux level at the flux max- 
imum and at longer wavelengths is matched by the model with 
T cff =175 000K. Models with T cS less than 170000K can be 
clearly ruled out, because the O v absorption edges become too 
strong. On the other hand, T e g higher than 180 000 K causes a 
strong flux excess in the region of the observed flux maximum. 

However, in one respect this model fails to fit the observa- 
tion. Its flux is distinctly too low at A < 100A. This could in 
principle be improved with a hotter model (200 000 K, along 
with a higher column density to fit the spectrum at long wave- 
lengths), but the flux exceeds the maximum by a factor of two 
in the range 100- 120 A (Fig.pl). This is the range where strong 
O vi line merging occurs and it is conceivable that, in case that 
we underestimate the O vi line broadening, stronger line wings 
can effectively block more flux in this region. Alternatively, 
yet unidentified opacity sources of other light metals may be 
responsible (e.g. Mg might be as abundant as Ne, and many 
Mg v and Mg vi lines are located in the EUV), but any progress 
must await data with better spectral resolution. 

Another argument that speaks for T e g slightly higher than 
175 000 K is the strength of the optical Ne vn 3644A line (the 
very strong EUV lines of Ne are much less temperature sen- 
sitive). A better fit is obtained at 200 000 K, although still, the 
equivalent width is too weak (Fig.H). This can be improved by 
reducing \ogg to 7.5 and/or increasing the Ne abundance (see 
below). 
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Fig. 3. The same models which were fitted to the EUVE spec- 
trum in Fig. are here compared to the ROSAT PSPC pulse 
height distribution. The 175 000 K model (dashed line) cannot 
fit the data but the 200 000 K model (full line) does, except 
for the softest channels (see text for discussion). Interstellar 
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respectively 
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As one can expect from the EUVE flux below 100A, the 
ROSAT PSPC pulse height distribution (120-30A) cannot be 
matched by the "cool" 175 000 K model (Fig. ph. A model with 
T c ff=200 000 K is needed to obtain a satisfactory fit (where we 
accept a model flux excess at the softest energies, which might 
be caused by the underestimated flux blocking in the 100- 120 A 
range as mentioned above). 

To summarize, we estimate T c ff= 170 000-200 000 K from 
the EUV and X-ray data, which compares reasonably well 
with the range obtained previously from the optical spectra 
(T ofr =150000-190000K). This indicates that H1504+65 ap- 
proaches or even exceeds the temperature of the hottest known 
PG1159 star (T cff =180000±20000K; RXJ0122.9-7521; 
Werner etal. 1996b). 

4.2. Neon abundance 

The EUV spectrum indicates a high neon abundance (2%), 
which is a factor of 20 higher than the solar value. We esti- 
mate a possible error of a factor of three, which is caused by 
the S/N of the data and the fact that the strong Ne vn lines are 
saturated. Our first detection of the Ne vn 3644A line in the 
Keck spectrum of H1504+65 (a previous attempt with the 3.5m 
Calar Alto telescope failed) indicates, that the actual abundance 
might be even higher. Fig.H shows several model fits to this 
line. Even under the assumption of a lower gravity (log g =7.5), 
which represents the lower limit of the optical analysis and 
which results in a stronger line, the relatively strong observed 
profile suggests an abundance of the order 5%. We successfully 
detected this line earlier in three other PG 1159 stars (Werner 
& Rauch 1994) and, as discussed in detail in that paper, the de- 
rived Ne abundance (2%) can be expected in matter processed 
by 3a burning. This high abundance can be rated as an indi- 
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Fig. 4. The Keck spectrum shows the Ne vn 3644A line and 
its strength confirms the high neon abundance deduced from 
the EUV lines. Left panel: Profiles from the two models fitted 
to EUVE and ROSAT data in Figs. | and |. With Ne=2% (20 
times solar value) the profiles are still too weak, however, they 
are also gravity dependent, as can be seen in the right panel 
(log g decreased from 8 to 7.5). It suggests that the Ne abun- 
dance might be even higher than 2% 

cation that the stars have lost their H- and He-rich envelopes, 
with H1504+65 being the most extreme case. 

We have already speculated that the He-deficiency in 
HI 504+65 points at the possibility that the star has burned 
carbon in previous evolutionary stages (Werner 1991). Hence, 
it might have been one of the "heavyweight" intermediate- 
mass stars (8Mq< M < IOMq) which form white dwarfs 
with electron-degenerate O-Ne-Mg cores resulting from car- 
bon burning. The extraordinarily high Ne abundance found in 
the present study and the unusually high mass of HI 504+65 
corroborate this speculation. In any case, the Mg and Na abun- 
dances can be as high as the Ne abundance (Iben etal. 1997), 
so that better EUV spectra, which can be provided by the fu- 
ture Chandra mission, should be used for further studies in this 
direction with a possible identification of other metal lines. 
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